here llie *ubsuipted Junctions are polynomials, lo apply llie general equation solving 
ehnique given at the beginning of section 2 to this problem, degree bounds must be 
dablished on the polynomials. 

The multivariate problem is even more difficult than the rational function problem, 
tere we are given a polynomial x 2 ,.. -* r ) and wish to know if there exist a 

(Mj, u 2 , * . •, u 5 ) and fi^x^ .v 2 ,. .., x r ),.. /t/xj, x 2 , - . x r ) such that 

4 /(xj, x 2 , - - -i x' r ) ~ 6l(/q(Xj, • . .» x r ),. . ., h s (x i, x 2 , .v r ). 

hough even partial solutions would be an aid to algebraic simplification and evaluation 
roblems, we know of no results in this direction. 

Others have attacked similar problems. Trager & Yun, (1976) presented an algorithm 
>r completing /ith powers of polynomials. That is, given a polynomial fix), they 
etermine if it can be written in the form a(g{x)) n + b for a and b elements of k. An earlier 
~‘rsion of the first decomposition algorithm given here was presented by Barton & Zippel 
! 976). 

Finally, we note some potential applications of our algorithm. In addition to solving 
olynomial equations it may be used for pre-conditioning polynomials which are to be 
:peatedly evaluated, e.g. in coding theory and numerical analysis. Also, the 
^composition of a polynomial gives the complete lattice structure of the subfields of the 
^nus 0 curve y=J[x). 

The presentation in this paper was substantially improved by comments from Gerry Roylance. 
upport for the preparation of this paper was provided under DARPA contract 
100014-80-C 0622. 
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We develop the theory for decreasing the depth of nesting in expressions that contain square 
roots. We show exactly when fourth roots enable denesting of expressions not denestable with 
square roots only. When we restrict our attention to denesting over the real numbers, we show 
in fact that no roots other than square roots or fourth roots are ever useful for denesting 
expressions containing square roots only, thus characterising the denestable expressions in this 
case. 

We then proceed to describe new algorithms that accomplish such denesting. 


Algebraic Simplification and Denesting 


Algebraic simplification is an important component of any symbolic manipulation system. 
In a general setting the simplification problem is undecidable (Richardson, 1968). For 
important special cases, however, simplification may be feasible. Here we are concerned 
with simplification of expressions containing nested radicals; in particular, we wish lo 
decrease the depth of nesting whenever possible. This process is called denesting. Two 
simple examples of formulas that denest are 


5 + 2^6 = 72 + ^ 

3 + 2^3 =iyi2(2 + yi2). 


Earlier work was done by Caviness & Faleman (1976) and Zippel (1977). Zippel 
showed how to denest certain expressions containing square roots using only square roots 
in the denesting, but left open the question of whether arbitrary roots are of any use in 
denesling such expressions. In section 1 wc develop the theory for denesting expressions 
that contain square roots. We show exactly when fourth roots enable denesting of 
expressions not denestable with square roots only. When we restrict our attention to 
denesting over the real numbers, we show in fact that no roots other than square roots or 
fourth roots are ever useful for denesting expressions containing square roots only, thus 
characterising the denestable expressions in this important case. 
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h< sections 2 6 we then proceed to describe new algorithms that denest these 
xpicssions. The algorithm of section 2 handles certain simple expressions of arbitrary 
.c:>ung depth, such as the following example of depth 3 


/ 16-2 n /29 + 2 v / 55-10 n /29 = 75 + ^11 -2^29. 


iection 4 adds the capability of fourth roots to the denested formulas. Section 5 extends 
he algorithm to handle rational combinations of such nested expressions, for instance 


/* +v / 3+’ 


6 + 373-710 + 6^/3, 


vhicli denests to 0, despite the fact that no individual term denests. Finally, section 6 
ucscuts an algorithm that denests more complicated doubly nested formulas, such as 


/(112 + -mjl) + (46+34 v /2) v /5 = (5+4^2)+ (3 + 72)75. 


li K is a field, then K(a { , . .., a„) denotes the smallest field containing K and a Jt .. a„* 
did is called an extension of K. If K(a lt . . a n ) ^ R, where R is the field of real numbers, 
hen tlie extension is called a real extension. The degree of K(a : , . . ., a n ) over K , denoted 
K{n lt . . ., a„):K~\, is the dimension of K(a u . . a n ) as a vector space over K, The reader 
s referred to Herstein (1975) for basic facts concerning extension fields. 

A formula over a field K and its depth of nesting over K are defined as follows: 


1. an element of K is a formula of depth 0 over K , t 

2. an arithmetic combination of formulas A and B is a formula whose depth over K is j 

max(depth(/4), depth(/I)), and j 

3. a root of a formula A is a formula whose depth over K is 1 +depth(/l). | 

lluoughout the paper the word root always refers to the principal (positive real) root. | 
V formula can be denested over K if it is equal to another formula of lesser depth over K. j 
f K is not explicitly mentioned, the words depth and denest refer respectively to depth ; 

)\ci <7, and dencst over Q, where Q is the field of rational numbers. * 


1. Basic Theorems for Denesting 

1 

In this section we ask when a formula of depth 2 over some field K can be rewritten as 
1 formula of depth 1 over K. Theorem 1, mest of which is implicit in Zippel (1977), shows 
hat if yfa + byfr can be denested using only square roots, then it can be denested 
ntioducing only one new square root. In addition, Theorem 1 gives a necessary and 
ulhcicnt condition describing when such a formula can be denested. t 

1 hcorem 2 is analogous to Theorem 1 and covers the case when fourth roots help in j 
leiU bling. Theorem 3 shows that roots other than square roots or fourth roots never help f 
fidi roots involved are real. * 


Ihiorem 1. Let ()c K and let < 1 , b,re K with yfr £ K. Then the following three statements 
ire equivalent. 

(1) y/a + by/re K(y/r, v 4i,. N /u fc ), for some a x 

(2) yfs(a + b^fr) e K(y/~r ), for some se K — { 0}. 

(3) 


a k e K . 


Proof. (1)=>(2): This part is by induction on k. 

Basis {k = 0): Choose s = 1. 

Induction (k > 0): Without loss of generality assume 


y/«ki K(y/r , yfa u .. Ja k . J. 


y/a + byfr = c+dyfa k9 

where c and d are in K(y/r, yfa[, ...» s/a k Z t ). Then 

a + by/r = (c 2 + d 2 a k ) -f 2 cdy[a k . 

Since a A-by/r, c 2 -\-d 2 a k and 2 cd are each in K(y/r, y/a/, .. y /a k S 1 ) but yfa k is not, 

= 0 . 

Case 1 (c = 0): Then y/a k {a + by/r ) = a k d is an element of K (y/r , yfa [,. . y/a k - x ). By 
the induction hypothesis there is an s in K — {0} such that y/sa k (a + by/r) is an element 
of K(y/r). Let 5 = sa k . Then y/s(a + by/r) is in K(y/r) and s is in K — (0}. 

Case 2 (d = 0): Then J a + by/r = ce K(y /r, y/a~ x ,. . y/a k - 1 )• By the induction 
hypothesis there is an s in K — { 0} such that yfs(a A- by/r) is in K^r). 

(2)-(3): 

Let 5(fl + ft7'") = (c + ^T*") 2 = (c 2 + d 2 r) + led Jr, where c, de K. Since ~Jr is not in K, 
sa = c 2 + d 2 r and sb = 2cd. If d = 0, then h = 0 and so y/a^—b 2 ? = ae K. Otherwise, 
eliminating c, 

s 2 b 2 J2 

sa = ^F +dr ’ 

or 

b 2 s 2 — 4ad 2 s + 4d 4 r = 0. 

Since seK, the discriminant \6a 2 d 4 — \6b 2 d*r must be a perfect square in K, or 


( 3 ) =>( 1 ): 

Let d = y/a 2 — b 2 r e K and s = 2(a-\-d). If s = 0, then a = —d = — yja 2 — b 2 r, so 
a 2 = a 2 — b 2 r, from which b = 0 follows, and (1) holds. Therefore, assume s ^ 0. Then 

7 a + by/r = 6 K(y/r, J s), 

27 s 

since 

/s + Iby/r'V 5 2 + 4 b 2 r + 4 bsy/r _ 4 (a + d) 2 + 4(a 2 — d 2 ) + 8 b(a 4- d)y/r 

V 2 y/s ) 4s 8 (a + d ) 

= \[(a + d) + (a — d) + 2by/?] = a + byfr. □ 

It is possible that a doubly nested expression y/a-\-by/r may not denest using only 
square roots but may denest if the fourth root of r is introduced. An example is 

74+372 = 7^(1+72). 




Applying the lest “ v a l -b 2 r e Q" of Theorem 1 to ^4 + 3^2 shows that it cannot be 
dcnested using square roots only. Theorem 2 gives the conditions under which the fourth 
root of r helps. 

TimoRHM 2. Let QcK and let a , b,reK with yfr $ K. Then the following three statements 
lire equivalent. 

(1) y/a + hyfreK(*/r, y/a[ . y/a k ), lor some a .. a k e K. 

(2) tf^ y fi+hj~reK( y /l) or y/iiy/a + by/r e K{y/r), for some s e K-{0}. 

(3) y/r(b 2 r—a 2 ) e K or s /a 1 -b 2 re K. 


Proof. (1) => (2): 

Assume y/a + b^/r is not in K(y/r, y/cii, . . \A*)i since otherwise y/s yJa-\- hyfr 

e K(y/r) follows from Theorem 1. Let a + by/r = (c 2 +cl 2 y /r) + 2cd */r, where c and d are 
in K(yfr, y /a l ,.. Furthermore, we can assume d ± 0 and 

f/rtK(y/r n ^/iii,...,y/ti k ), since Ja + byfr is in K{^Tr %s fa[. .\Ac)~ 

K(y/r, y/ci^.. yja k ). Hence c = 0 and ,/« + fc,/' 7 = dtfr. Multiplication by tfr yields 



/br + dy/r. 


This means y/br + (iy/r is in K(y/r, v /^7,. • *, y^) : By Theorem 1, y/s(br -f is i n 
K(y/~r) for some s in K — {0}. That is fjr yfsyfa + by/r is in K(y/t m ). 


(2) =► (3): 

If I /i’/r is in Kf^/r), then Ja 1 -h 1 r is in by Theorem 1. Otherwise 

*/r v /.v n /« + /) v 4 is in K(y/r), or equivalently s fs(br + ajr) is in K(yfr). ByTheorem 1. 
y/b^r 2 — a 2 r is in K. 


(3) =>• (1): _ 

If yfa 2 -b 2 r is in K , then ^/u + b^ is_ in K( y/>\ y/s), by Theorem 1. By Theorem I. 

in K implies Jbr+Oyfr = tfrja + byfr is in K(y/r, y/s) for some s in K, so 
y/a + by/r is in K(f/r, y/s). □ 

Theorems 1 and 2 will be used in algorithms that denest formulas. Theorem 2 raises the 
possibility that roots other than square roots may be useful. In Theorem 3 we show that, 
with the exception of the fourth root of r, no additional roots are useful and that 
Theorems 1 and 2 cover all possible denestings. 

Whereas Theorems 1 and 2 hold for arbitrary extensions. Theorem 3 holds only for real 
extensions of K. For example, yf\ + yj— I docs not denest using only square roots and 
fourth roots, but can be denested using an eighth root: = 7-4. One may 

question whether this is really a simplification, or whether one would prefer to treat *f-4 
as \J\f\f ~4 and denest to v /l + N /—"l. 

To simplify the proof of Theorem 3 we first prove five lemmas. Given a real extension 
of K in a particular normal form, Lemma I states that for any ,v in K with y/x in the 
extension, y/x can be expressed in a simple form. Again, given the normal form. Lemma 
2 states that certain roots of elements of K are not in the extension. Lemmas 3 and 4 


establish that ccilatn reul extensions of K can tie placed in die noiiiuii tonu. i '-"""“ - 
deals with the degree of an extension obtained by adding a prime root. 

For similar results in the case when the underlying field K is Q, the reader is referred to 
Besicovitch (1940) and Richards (1974). 


Lemma I . Let K be a real extension of Q , and let a„ ..., a n+ , be positive elements ofK. Let 
. r„ 5 2 he powers of 2, where r, =/= 1 for 1 < i ^ n. Suppose that 

[K(7b.^):K] = fl r, 

i = 1 

and that the r { are minimal powers of 2 such that y/a n +7 hs Then 

r { = r 2 = • • • = r n = 2 and y/a„ +l = a y /a l • • ■ 

for some a in K. 


Proof. The lemma is proved by induction on n . 

Basis (n = 0): Let a = Ja n+ 

Induction (n> 0): Let L = K(a\' r \ ..., fl 2/r ") and write y/a/7, =ba l f r " + c, with 

b and c in L. By the minimality of r„, b # 0. By squaring, 

a n+ , = b 2 a 2lr " + c 2 + 2bcn„ I/r ". 

Since b^0 and a}/ r " is not in L, we conclude that c = 0. Assume now that r„ = 2. Then 
\ tt n+l a„ = ba„, so 

r . a »-r')- 

Therefore, by the induction hypothesis, each r,= 2, and there is a e K such that 


'a„ + ,a n = ay/a, 


>7 + i = tia;'y/a, 

Thus assume r„ > 4. Let L= K(a \ lr ',..., flJ'.'T ’, a*J r ") and write 

= bn,! /r " = (d+en„ 2/r "K /r ", 

where d, eel. Again, by squaring 

a n+ , = (d 2 + e^)aff"+2deaT". 

Since [L: L] = 2, a;f" is not in L, and therefore d 2 + e 2 a*f'" = 0. This contradicts the 
hypothesis that a„ is positive, so the assumption r„ > 4 is false. □ 

Observe that the hypothesis that K is a real extension is needed, since 

72 = (1 -7^)7^ e T). 

Lemma 2. Let K be a real extension of Q and let a be a positive element ofK. Let r >2 be a 
power of 2. Suppose [fC(n ,/r ): K] = r. Then n ,/2r 4 K<a‘"). 


Proof. Suppose that a' l2r is in K(a l ' r ). Write a' 12 ' = bn 1/r + c, where b and c are in K(« 2/r ). 





juaring \\e get 


a {,r = b 2 a 1,r -\- c 2 + 2bca l,r . 

Since a i,r is not in K(cr' r ) we conclude that 2bc~\ k so b 2 a 2lr + c 2 = 0. Thus either 
b = c — 0 or a 1,r = — (c/b) 2 , in either case a contradiction. □ 

1 i mma 3. Let K be a real extension of Q. Let n lf . . n k be powers of 2 and let 

k 

« = n 

i= 1 

Let tij, . . ., a k be positive elements of K. Let 

l = k{^ u .. 

If [L.K] n, then we can find positive b { , . . b k in K , and integers m l9 .. m k , each a 
power of 2, such that 

L = KCjb ~„..^//Tj and m = f{ m, < n. 


Pkuof. Renumber so that Mj ^ n 2 ^ ^ n k . Consider the sequence of extensions 

K, K{^), K(*/a 7),..K(T^7), K(-v^7, 7^), K("/^, tfT 2 ) . 

C iMibider the first step at which the extension is not proper. Then by Lemma 2 we must be 
audiiig yja { for some i ami 

yfeii e K(a\ ,n \ . . ., a}L "[~'). 

by Lemma 1, y /a i = b > /a il • • ■ a ir for some be K and i u ..., i r e {1,2,.,., i- 1}. Let 
L = K(a\' n \ .. a\L n \~\ b 2/ni , a\^\ .. al /nk ). 


a i\ ’ ^i, 


it is clear that Zv = L and m = a/2. □ 


Li mma 4. Given L as in Lemma 3 wc can find positive b l ,...,b k in K and integers 
m i, . . each a power of 2, such that j 

L= KCyr i ,...,’*/F k ) and [L:K] = f\ m, j; 


Proof. This follows from repeated application of Lemma 3. □ 

Lhmma 5. (Abel (Schinzel, 1982, p. 91)). Let L be a real extension of the rationals , p an odd 
prime , b a real number , <iw/ b p an element of L. Then, provided b<fL, [L(b); L] = p. 

Proof. Any irreducible factor of x p -/? p over the reals is x-b or of the form 
\ 2 — {(o + fo)bx + b 2 , where to and Co are complex conjugates, (This is because the complex 
roots of x p -b p = 0 are of the form cob, where co is a pth root of 1, and because 

(x — cob)(x — a?b) = x 2 — (cu -f Co)b + b 2 ). 


Suppose x p -b p could be factored over L. Since L is a real extension, any proper factor of 
x p -b p over L must be a product of some proper subset of its factors over the reals, rnd 
hence have a constant term of the form b\ 1 ^ i < />. Since p is prime, b ** e L implies b -- L, 
a contradiction. Thus we conclude that x p -b p is irreducible over L and horce 
[L(b ): L] - p (see Herstein, 1975). □ 

As an example of the importance of the real assumptions, let 


1 +J — 


and p = 3. Then b p = - 1, but Q(b) = 3), so [Q(b) :Q] = 2. 

Theorem 3. Let K be a real extension of Q and let a,b,reK with yfr K . Let 
«i,. ■ n k ^ 1 and a l9 ..., a k e K be positive. If 

y/a + by/r e K(ffa[, ■ ■ ■, l/th) 

then _ _ _ 

y/a + by/reK(*/r, Jp ~ u ..yfp k ) 

for some positive p l9 ..., p k in K. 


Proof. Let n u . . n k be minimal such that y/a + by/r is in K(a\ ,n \ .. ., a\ lnk ). Suppose n k 
were divisible by an odd prime p. Let L = K(a\ ln ‘,..af "\~By the minimality of n k , 


L(ar k ) $ L(af"\y/a + l 


; L(af'). 


But yfa + byfr ): L«"-)] is 2 or 4 and [L(a‘ /n -): L(up")] = p (by Lemma 5), a 

contradiction (see Herstein, 1975). Thus each n t is a power of 2. 

Now let L = ZC(a} /M1 ,.. ., af nk ). By Lemma 4 there exist positive b { . b k in K and 

integers m l ,...,m k , each a power of 2 and each at least 2, such that 
L = K(b\ ,m \ . . ., b k lntk ) and 

[L:K] = n m , 

/= i 

Assume m l ,... 9 m k are the least powers of 2 such that yja + bj~r is still in 
K{b\‘ m \ ..., bl ,ntk ); if not, de crease th em accordingly and redefine L. Since ^/reL, 
Lemma 1 states that J~r = c y /b l •" • bj for some ceK. Assume that m i ^ m 2 ^ ^ m p 

possibly renumbering the b { . Let r i =b i ---b J . Now L = K(r\ ,m \ b l 2 lm \ . . b k lntk ). 
Suppose by way of contradiction that m { ^ 8, or i > 1 and m { ^ 4, and denote r} /mi (in the 
first case) or b\ lm ' (in the second) by p 1/m . Let L=K{b l f m \ . . ., b l k /mk ) (in the first case) and 

L = K(r\> m \ . . b\LT\ b}i m ‘ +i ,. . ., bi> mk ) 


fin the second). Observe that y/r = c y Jr l is in L(p 4,m ). 

By the minimality of m, yja Lb Jr is in L(p l/m )- L(p 2,m ). That is, 

yfa + by/r = c + dp 1,m , 

where c and d are in L{p 2lm ) and d #0. By squaring, 

a + by/r = c 1 + d 2 p 2,m + 2cdp l,m . 



Since all bul tin. last Icrm irr//) 1 arc in L(p ~"'), it follows that cd — 0, so c — 0. Wrile 
Ja + h^fr = dp'"" = (e + fp 2m )p' ,m , 
where e and/are in Up 4 '"'). Again, by squaring, 

a + l w r r = ( e * +/V"V / " + 2c/p 4 "". 

All subformnlas of this equation are in /„(/> 4 "") except for p 2/ '\ since 

IL:K] = fl »>,• 

I = 1 

Thus e 2 +/ 2 /j 4/ " -0 implying that e =/= 0 or = -r//' 2 , a contradiction. Tims, is 
in {2, 4} and m, = 2 for i > I, i.e. 

1 1 \fb~2' • ■ y/bk)- 

=y;- 4'V,, y^yy e y?, y^;. vm- □ 

2. Denesting Using Square Roots Only 

The goal of the remaining sections is to present new algorithms that, whenever possible, 
rewrite a given formula with positive radicands as one that has a lower nesting depth over 
Q. (The algorithms will sometimes, but not always, denest formulas that have negative 
radicands.) Consider the following tower of fields: 

K ( > = Q and, for all i ^ 0, K H , = K^y/r^ 

where r t eK livery element of K t has depth at most ? over Q. Consider y/r a , where 

r n e K n . This section presents an algorithm, based on Theorem 1, that rewrites JV n as a 
formula of depth at most n over Q involving only square roots, whenever such a formula 
exists (even if the formulas have negative radicands). 

Although the fundamental ideas appear in this algorithm, it is incomplete in three 
respects. Firstly, it does not permit the fourth roots of Theorem 2; this is rectified in 
Section 4. Secondly, the algorithm denests only single radicals from K n + i rather than 
arbitrary elements of K fJ , this is rectified in section 5. Finally, not every formula of 

depth n over Q can be expressed in the form of elements of /C H ; in particular, elements of 

K n can be written using only n distinct radicals. Section 6 rectifies this only partially, by 
describing how to denest arbitrary formulas of depth 2 over Q, and explaining why this 
technique does not generalise to depths greater than 2. 

Let r„ = u + b~Jr n _,, where a, b , r„_, e K„_,. There are two ways in which J7,, might 
denest using square roots only. The first way is that v n j might denest. The second 
possibility is to find some field K containing only elements of depth n — 1, containing tu 
and r M _ l5 and in which a 2 — b 2 r n -i is a perfect square. Then y/r~ n denests by Theorem 1. 
In order to find the K in which ^ 2 -/7 2 r /I _, is a perfect square we attempt to denest 
J a 2 - b 2 r n 7^ Tf we are successful we will find a field K, all of w hose elements have depth 
at most n — 2 and an element s in K such that y /a 2 — b 2 r „-j is in K(y/s, >/ r «- 2 )- The 
desired K is then K(y/s, N /r„- 2 ). 

Since y /r n . l and y/a 2 -h 2 r~i each have depth less than y/r„ this suggests a recursive 
algorithm: to denest y/r, n try to denest each of yfr^l and v /<i 2 -f> 2 r n _ 1 . If either 


succeeds, then yjr n can be denested, If neither succeeds, y/r n cannot be denested using 
square roots alone. 

This algorithm has a serious shortcoming: it requires exponentially many recursive 
invocations. However, closer inspection reveals that there are not 2 1 distinct arguments at 
the ith level of recursion, but only i+1, as recurs repeatedly. Continuing the 

previous discussion of denesting ^/r^, the first level of recursion attempted to denest 
yfr^ and Jic - b 1 ^^ Each of these radicands is in K n . u so 

= e + djr^. 2 , and a 2 -/> 2 r„_, = e+fjr^, 

for some c,il,c,/,r,. ; eK„.j. Thus, the only formulas t hat need b e denested at the 
second level of recursion are ,//•„_ 2 , sjc z — il 2 r„_ 2 , and c 2 —/ 2 r„_ 2 - In general, only 
one more argument will be added at each level of recursion. 

This motivates the algorithm of this section. We assume for the remainder of this paper 
that there is an appropriate data type formula for representing formulas. The recursive 
algorithm takes as input 2 integers, n and m, and a list of m formulas, each of depth n over 
Q. The algorithm attempts to denest each of these, returning for each input formula either 
a denested version, or the input formula itself if it was unsuccessful. 

Notice in the algorithm that it is assumed that each inp ut formul a shares the same 
value of y/r if it ^ 2; that is, the ith formula is of the form y/of + bfy/r. 


function denest ( integer , n, m, 

array [1 ... m] of formula nested) 
returns array [1 ... m] of formula ; 
begin 

declare N, d: array [1 ... m+1] of formula ; 
declare s: formula ; 

If n = 1 

then for i from 1 to m do 
</nested[i] = y/x 2 for some positive xeQ 

then denestJj] <- x 
else DENEST[i] <- nested[i] 


else begin 

assume nested[i] = yj a t + b iN /r, 1 < i ^ m; 
for i from 1 to m do N[i] <—^/aj 2 — b 2 r; 

N[m + l]4- N /r; 
d<-DENEST(n — 1, m-fl, N); 

//N[m+1] ^dfm + l] comment : y/r denested; 
then for i from 1 to m do 

denest[Q 4- yai + bj'dfm + l] 
else for i from 1 to m do 
if N[i]=d[i] 

then DENESTfi] <- nestedfi] 
else begin 

s 4 -2(a i +d[i]); 


DENESTfj] 


s+2b iN /r 

2 y/s 


end 

end DENEST. 


end 


Consider the following example due to Shanks, see Zippel (1977), the object of which is to T 
denest | 

^16-2729 + 2^55-10^29 f 

Here I 

,nH Kl = Qi ^ 9) I 

K 2 = K,(V55—IO 729 ), I 

satisfying 16 — 2729 e K { , 2 e K t , and 55— 10729 e K t . The first call is I 

denes r(3, 1, [716-2729 + 2755-10729]). I 

This produces the two formulas > 

7(16-2729) 2 -2 3 (55-107291 = 7152-24729 and 755-10729, 
resulting in the second call of 

df.nest(2, 2, [7152-24729, 755- 1°729])- r 

This in turn produces the three formulas k 

v /l52 2 -24 2 x 29 = ^6400, y/55 2 - I0 2 x 29 = 7 T 25 , and ^29, 

resulting in the third call of denest(1,3, [^6400, 7T25, 729]). The basis portion of the j" 
algorithm succeeds in denesting the first of these, returning the list [80, ^/29]. The jj 

invoking procedure can use this to denest its first input yj\ 52 — 24^29 as follows \ 


s = 2(152 + 80) = 464, 
7^ 464-2x24^/29 


= 2729 - 6 . 


(The last simplification would not be done by the algorithm, but it makes the example 

clearer. Notice that it does not afTect the depth of nesting.)_ 

Thus, this invocation returns the list [2y/^~6, y/55^-Jo^/S]. Finally, the initial 
invocation completes as follows 

.V = ((16—2729)+ (2729 — 6)) = 20, 

7 ^ 729+2755 — 10729 = 20+2 x 2 77 ^ = 75 +7* 1 - 2729 . 


returning the singleton list [75 + 711—2729]- 

The correctness of this section’s algorithm is the subject of the next theorem. 

Theorem 4. Let m^\ and h>1. Given m formulas 7/7 yJTi > • • •» \[Tm-> where 
fu /21 • • •» /m e K„_ 1 > procedure denest will denest y Jf i correctly , whenever y Jf i can be 
denested using square roots only, for all 1 ^ i ^ m. 

Proof. The proof is by induction on n. 

Basis (n = 1): f e Q , so 7 denests if and only if f is a perfect square. 


* 



I 


Induction (n > 1): Let f = a { + b iy /r n Z 2 , where a h b h r„_ 2 e A7_ 2 , and suppose JJ- X can 
be denested to depth n — 1 over Q using square roots only, that is. 


/ a i T biy/r n _ 2 £ K(77, 7^2» * - *> 7 / ^*)» 


for some K all of whose elements are of depth at most n — 2 over Q, such depth achieved 
using only square ro ots, a nd some c,, c 2 , ...,c k eK. We can assume that /C„_ 2 ^ K. By 
Theorem 1, either 7 r »-2 6 K or yjaf — bfr n _ 2 e K, that is, one of these two formulas 
denests. Each of these two radicands is in K n _ 2 , and each of these two formulas is passed 
in the recursive invocation so, by the induction hypothesis, denest will denest one or both 
of them corre ctly. If yjr n _ 2 is denested to reK , then surely y/aq + bj is a correct 

denesting of 7 <+ +Otherwise, suppose J a 2 -b}r n _ 2 is denested to e K. Then 
the correctness of the denesting 


7 + ^7 r "-2 


h — 2 1 yfi). 


where s = 2(<7 t - + d,.) is supplied in the proof of Theorem 1, part (3)=>(1). □ 


What remains is a discussion of the running time of the procedure denest. Inspection of 
the procedure reveals that, as long as the formulas have size polynomial in n, the number 
of arithmetic operations will also be polynomial in n. Unfortunately, in general the 
formulas have size exponential in n , as is clear from the last two lines of the procedure: 5 
depends on d\f_ J, the result of the recursive invocation, and there are two occurrences of s 
id DENESTfi], the result of the current invocation. We alleviate this problem in an 
(unfortunately) asymmetric way: although we require the inputs to be represented as 
trees, we represent the outputs as acyclic directed graphs (i.e. straight-line programs). 
This is unfortunate in the sense that the output cannot be given as an input for another 
level of denesting. Notice also that, unlike the arithmetic complexity, the bit complexity is 
not polynomially bounded, as the number of bits doubles at each recursive call due to the 
squaring involved in yja 2 — b 2 r. 


3. Finding a Perfect Square Among Products of a Set of Numbers 

Before proceeding with the next algorithm we present a straightforward technique for 
finding a perfect square among all possible products of a set of rational numbers. The 
result is used as the basis for the recursive algorithm in the next section. 

Let r u .. ., r„ be rational numbers. Let r t - = r^f, where d { is the denominator of r f . 
(Note that r f is the product of the numerator and denominator of r,.) Then Ylies r f is a 
perfect square if and only if r t is a perfect square. Thus, without loss of generality, let 
r,. r n be integers. 

In order to find a perfect square among products of r u r 2 , . . ., r n the greatest common 
divisor of each pair is determined, resulting in a set of integers pj that are pairwise 
relatively prime and for which each r t -, 1 ^ i ^ n, can be written as a product of pjS. 
Remove any pj that is a perfect square. 

Next, for i = 1, 2.w, replace r ( by p { a h where a { is a product of pf s, j > i. (This may 

involve renumbering the p- s.) In order to do this, if r ( contains p jt j < i, then for the least 
such j replace r, by r,^, remove any p k from this product, and iterate the process until r i 
contains pj only for J>/. (At this point renumber one such pj to be p h if p t does not 



already occur in r,.) Il is straightforward to verify that such a replacement does no! 
change the existence of a subset whose product is a perfect square, and that such a subset 
exists if and only if some r, eventually takes on the value 1. The subset can be identified 
easily by retracing the steps. 

The following application of this algorithm will be used in section 6. In the case where 
all r,-> 0 we can determine if y/r„ is in Q{ s fr [ t .. .,^/r /l _ 1 ) by locating perfect square* 
among products of r if . . r„. On encountering a perfect square, if r„ is involved, then 

is in Q{yfr u ..., v /r^_ j). Otherwise delete some r, from the product, and repeat the 
process until y fr u is discovered to be in Q(y/r^, .. yjr n _,) or there are no perfecl 
squares remaining. In the latter case s fr n is not in Q(yfr [,. . y/r n -i) by Lemma 1. 

4 . Dencsting Using Arbitrary Roots 

Consider again the tower of fields presented in section 2 (but in this section restricted to 
be subfields of the reals) 

Ko = Q, 

and 

for all i ^ 0, K i + , = 

where e is positive. The goal of this section is to generalise the algorithm of section 2 
so that it rewrites its input y/r~ n as a formula of depth at most n over Q , whenever such a 
denesting is possible, yet still has arithmetic complexity polynomial in n. By Theorem 3. 
the only useful additions to the algorithm of section 2 are the fourth roots of Theorem 2 
(provided all subexpressions that arise arc real). 

Let r„ = a + /\/r„ _ { , where a, b, r„_, e K n ~By Theorems 1-3, yfr n denests if and 
only if 

(i) denests, or 

(ii) y/a 2 ~b 2 r„- l denests, or 

(((•) denests. 

Unfortunately, the obvious generalisation of the algorithm of section 2 has exponential 
complexity, since each formula y/a -f b y jV n _ j in the input list gives rise to two formulas. 
yja 2 —/) 2 r n _, and yf—r lt fa 2 — b 2 r n _ j) in the recursive input list, in addition to the one 
extra formula y/r„_ i . 

The solution is to deal with the formulas of the form yj— r n _fa 2 — b 2 r n - { ) implicitly, 
rather than generating them explicitly, by observing that the radicand can be expressed as 
the product (with a possible sign change) of the radicands of yjr n _ j and y/a 2 — b 2 r n _ x . 
Thus, the algorithm of this section will invoke itself with exactly the same lists of formulas 
as the one of section 2, but rather than trying to denest each individually, it will determine 
whether there is a subset of the formulas whose product (with a possible sign change) 
denests. The justification for looking at all subsets is given below in Lemma 6, which 
states that a product 

k 


i = 1 

of input formulas denests if and only if yjr denests, or the corresponding product 

k 



ol formulas denests, or 


denests. 

In the basis, when all formulas have depth 1 over 0, the algorithm applies the 
procedure of section 3 to determine whether any subset of its input radicands has a 
product (with a possible sign change) that is a perfect square. If it finds such a subset it 
will identify it in the global array factor by setting 


„ _ fl, if the tth formula is in the subset 
factor[Q = < 

(0, otherwise 

The pending invocations will then use Theorems 1 and 2 and Lemma 6 to denest the 
product specified in the array factor. 


function denest (integer n, m, 

array [1 . .. m] of formula nested) 
returns formula ; 
begin 

declare N: array [1 . .. m+ 1] of m formula; 
declare d, a, b, s: formula ; 
f/n= 1 

then Use the procedure of section 3 to find a perfect square x 2 among products of the 
m input radicands. Before entering the radicand of nestedfm] into any product 
with at least one other radicand, it should be multiplied by -I. Record the 
formulas entering into the product by placing l’s into the appropriate cells of the 
array factor. Return the formula x. If no such perfect square is found, set factor 
[i] = 0 for all i and return the formula nested[mj. 
else begin 

assume nested[i] = b i% /r, 1 < i < m ; 
for i from 1 to m do N[i] <- yjaf — bfr; 

N[m+ 1 ]«— -v/r; 
d <- DENEST(n — 1, m 4-1, N); 

if Tactorfi] = 0 for all 1 ^ i ^ m comment : only could have denested; 

then denest <- y/a. m + b m d 
else begin 

, - -= «n 

define a, b by ^/a + b^/r = factor[i]* nested[i]; 

if factor[m] = 1 and factorfi] = 1 for some 1 < i < m 
then begin 
a<—a; 
b<—b 
end; 

comment : y/a + by/r denests; 
if factor [m+ 1] = 0 
then begin 

comment : apply Theorem 1; 
s<-2(a + d); 


DENEST <— 
end 


s + 2b v /r 
2,/s 


else begin 

comment: apply Theorem 2; 
s<-2(br + d); 


DI-NEST 


end 

end 

end DiiNiiST. 

Consider denesting the formula yfi 4 -^/s 7 . The depth of this formula over Q is 3, 

so the global declaration needed is 

declare factor: array [1.3] of { 0, 1}, 

and the first invocation is denest(3, 1 , [yj 3 + ^5 + 2^/7]). This gives rise to the two 
formulas 

s/tf — (5 + TJT) = and 75 + 2^7, 

resulting in the second invocation r>ENEST(2, 2, [^/4 — 2^/7, + 2^/7]). This in turn 

produces the three formulas 

v / 4 2 -2 2 x'7 = y^T2, V5 2 -2 2 x7 = y^3, and ^7, 

resulting in the third invocation denest(1, 3, [y/— 12, ^ — 3, T?])- The algorithm of 
section 3 determines that (— 12)( —3) is a perfect square, sets factor to [1,1,0], and 
returns 6. In the second invocation, a = 8 and b — 2, since 

>/(- 1X4-2/7X5 + 277) = V8 + 2 V / 7. 

Since factor[3] = 0 (i.e. y/—r was not involved in the product), s = 2(8+ 6) = 28 and 

so 1 + y/l is returned to the first invocation. In this invocation, a = 3 and b= 1. Since 
factor[2] #0 (i.e. J— r was involved in the product), a fourth root is needed as follows 

s = 2((5 + v /7) + (l+ > /7))= 12+6/7, 

+ 7^ + 2 = ((12 + 6/7) + 6/5 + 2/7)/(2/12 + 6/7/5 + 2/7) 

_ 6 + 3/7 + 3 /5 +2/7 

/l 2+6/7/5 + 2 / 7 ' ; 

The correctness of this section’s algorithm hinges largely on the following lemma. 

Lemma 6. Let k ^ 1 and n ^ 0. Let r e and, for all 1 ^ i ^ A:, let a h b t e K„. Then j 

* J— - r 

11 v T bjyjr 

i- i 

denests using only reals if and only if at least one of the following holds: j 

(0 Vr denests using only reals , or 


(//) P[ ^/a, 2 — bfr denests using only reals , or 


(w) yf’-r PJ y/af^bfr denests using only reals. 

i= 1 

Proof. Let a, b e K lt satisfy 

fl + fc/r = n (<J, + /j ; /r). 

i -1 

By Theorems 1-3, ^/o + ^/r denests using only reals if and only if at least one of the 
following holds: 

(i) y/r denests using only reals, or 

(ii) Ja 2 — b 2 r denests using only reals, or 

(iii) \f—r yja 2 — b 2 r denests using only reals. 

Hence it suffices to show that, under the hypothesis y/r $ K „, 

a 2 — b 2 r = P[ (af-bfr). 

i= 1 

This is done by a straightforward induction on k. □ 

The correctness of the algorithm follows by taking m = 1 in the next theorem. 


Theorem 5. Let m ^ 1 and 1. Given m formulas y / j \, y/f^, .. yjfm, where 
fofzi - • -ifm £ K n - i> procedure denest will correctly denest one of the following products in 
some real extension of Q , whenever possible: 

(0 s/Jm . or 

(ii) n y/f" f or some nonempty I £ {1,2,..m- 1}, or 

iel 

(ih) J~f m n y/f I* for some nonempty I c {1, 2,., m - 1}. 


Proof. The proof is by induction on n. 

Basis, (n = 1): This is precisely what the algorithm’s basis does. 

Induction (n > 1): Let one of the products that denests be rewritten as 

yja T byjr n ~ 2 — PJ \J c { + d [y J r n _ 2 , 

iel 

where r n _ 2 e K„_ 2 , and </,* e /C„_ 2 for all j e/. By Lemma 6, this product denests using 
only reals if and only if at least one of the following products denests using only reals: 

(i) y/r//>, or 

(>') n //-<//„-2. °r (4.1) 

iel 

(iii) / —r„_ 2 FI \/ c ?- d f r n-2- 


Since N /r fl _ 2 and y/cf -dfr n _ 2 (for all ie/) are among the formulas passed in the 
recursive invocation, the induction hypothesis says some such product will be denested 
correctly in some real extension of Q\ assume it is one of the three products in (4.1). If 

7 




</ = y/r„- 2 denes ted and is real, then surely 
is a correct denesling and is real. If 

1 1 = n 7‘F = ^77 = Vn (ct-dir,- 2 ) = v4T 2 -* 2 r,_i 

ie l iel 

denested and is real, then the correctness of the denesting 


V"+V r .-2-— 2 ^r. * 

where s = 2(«+d) is as in Theorem 4. Notice that ^/s must be real, since all other 
subexpressions are real. Finally, if 


s + 2byjv, t - 2 


(I = sf~r«- 


denested, then the correctness of 




/ “T~ 5 + lity r n-2 

'a+byJr „- 2 — r a/ - 

2Jsyr n - 2 


where s = 2{br n _ 2 + cl) is most easily seen by squaring both sides. That y/s is real is 
proved as before. □ 


5. Denesting Rational Combinations of Radicals 

So far we have developed the theory needed to denest certain single (possibly multiply 
nested) radicals. Suppose we have a linear combination of several radicals, none of which 
denests individually. Is it possi ble th at t he entire expressio n could denest? The answer is 
yes. For example, none of yf \ + N /3, >/3 + 3^/3, and yf\ Q+6^/3 denests, but 



equals 0 and hence clearly denests. 

Let K be a real extension of (?, and let /,*, fl,*, b i9 € K for !</</?, with 
Cl . + b is /r] > 0. (In fact, this last condition turns out to be unnecessary: although we need 
it in order to apply Lemma 1 in what follows, a weaker version of Lemma 1 involving 
only square roots could have been used.) 

Suppose that the linear combination 

ft 

E ' 

»= l 

denests in K. Then it can be denested as follows. First denest any individual radical that 
denests. As we shall see, it suffices to determine, for each pair of remaining radicals 
Kfii and \/*/ + fy>Z7’^eir product d enests in K. If the procedures of sections 
2-4 show that + biyfri y/ (l j +denests as keK, then replace 

hy/ a i + fy\A« + + fy*>/0 ^ 






and iterate the process of looking for a pair of radicals that denests. 

If at some point no product of a pair of radicals denests, then we claim that the entire 
linear combination could not denest. For let I = {/j, i 2 , • • •> 4,} a maximal set that 
satisfies 

[K(7 fl i.+*i.7 r «.> • ■ - 7°u + W r /J : *3 = r • 

By Lemma 1, for every j in {1, 2,. . ., /j} — /, 


y/ a j+ b jyfrj = n II V + 

ielj 

for some / y c / and some n e K. But the 2 m distinct products are linearly independent, so 
the linear dependence can occur only if I } - I k for some j^k, in which case 
yj + \/ a k + b k \fi\ denests in K. .... 

Note that linear combinations of radicals is as general as rational combinations, since a 
product of radicals is itself a radical, and quotients of radicals can be eliminated using 
cor>gates. 

Consider the example 

V / 1+V3 + V 3 + 3 V / 3-V /i0 + 6 n/ 3 - 

The product 7 1 +73^3 + 3^3 denests as TH^T 3 )- Thus we simplify the 
expression to 

^1+73(1 + % /3)- v / 1° + 6V3. 

Next the product J~\ +73^10 + 6^/3 denests as 4+ 2^/3. Thus we can simplify further. 

71 ^ + 73 + 373-710 + 673 =^ 1+73 


— -v/l 4-%/3 


. ,TVJ i+V 3 J 

(4 + 273) — (4+273) ' 

l+T 3 


6. Denesting Arbitrary Formulas of Depth 2 

In this section we consider a slightly different tower of fields 

Qo — Qi 

and — 

for i ^ 0, Q 1 + 1 = 6<(v r »)’ 

where r t e Q is positive. Using the procedure of section 3, we can assume that y/r t 4 Qi for 
all i. What makes the field Q h different from the field K h considered in section 4 is that any 
element of Q h has depth at most 1 over Q , for all h. 

Let f e Q h . In this section we describe an algorithm that denests y/f using square roots 
only, whenever such a denesting is possible. As a formula over Q, / is the linear 
combination of as many as 2 h different square roots, so the algorithm of section 2 is 
inapplicable to Jj\ As a formula over Q h -i, however,/can be written as 

/= a + \Aft-i, 

where a, be Oh-,. Thus, we can use the algorithm of section 2 to denest 


where a y bGQ h - l . 



v'7’- \/a + by/r h - ! over Q h _ lt (The reader is referred back to page 170 for a reminder of 
tlie definitions of “depth over Q h -f' and “ denest o ver Q h _ 

Since , is not an element of ^/a + bJnii has depth 2 over O. When the 


Since y/r h . r is not an element of Qh-r._\6 , + ^v r *-J has depth 2 over When the 
algorithm of section 2 denests y/o + by/r/- i over Q h _ 1? it will be rewritten as a formula 
of depth 1 over (?*_,, that is 


V 'H- byjr k .., e Q h _ , (y/r h _ j , N /a ,, V a*), 

for some a u a 2 ,. .e Q h - t . Notice that this “denested” version, like yfii + by/r h _ ,, 
still has depth 2 over Q, since a t , a 2 , . . a k in general have depth 1 over Q. The progress 
we have made is that we have eliminated one radical, y/r h _ ln from all depth 2 
subformulas, since a lt a 2 , . . a k are elements of Q/, — i, but y/r h _ j is not. Thus, we can 
denest yfo/ 9 yfci 29 * ■ yj'th recursively. 

Given a positive fe Q h , the procedure that denests y/f over Q will return deQ h and 
positive se Q satisfying y/f— dy/s. The procedure itself is as follows: 

1. If /i = 0, then return d if f=d 2 for some d e Q , and fail otherwise. If h >0, then write 
f=a + byfr h - u where a.beQ^^. 

2. Compute a 2 — b 2 r h _, e Q h -i. 

3. Recursively denest y Ja 2 -b 2 r h _ l over Q. If this fails to denest, then y/f fails to 
denest. Otherwise, the recrusive invocation returns eeQ h ^ i and seQ satisfying 
y /a 1 -h 2 r h _ i =e y /s. 

4. Use the procedure of section 3 to determine whether y/s e Q h - Y . If not, y/f fails to 
denest. If so, the procedure expresses y/s as an element of Q h - V 

5. Let d = e^Ts eQ h ^ i , and t = 2(a -t- d) e Q h _ j. Then 


If t e Q, then return 


and halt. 


1 + 'bill 


e Q h and t , 


6. Recursively denest J~t over Q. If this fails to denest, then y/f fails to denest. 
Otherwise, the recursive invocation returns jeQ h _ l and u e Q satisfying /I = jyfu. 

_ /- r~r - t + 2by/i\l 7 f + 2/\/r7_ 1 , _ 

7. y/a + by/r h _ !=— - Y ^ —, so return --- e Q h and // e Q. 

2jy/u 2uj 


e Q h and u e Q. 


Unlike all of the previous procedures, this one has exponential arithmetic complexity 
(measured as a function of h). Notice, though, that some formulas/in Q h may have size 
2 h ; for these, the arithmetic complexity is polynomial in the length of/. 

The correctness of the procedure is the subject of the next theorem. 

Thi-orfm 6. Given a positive fe Q h , the procedure will denest y/j over Q correctly whenever 
^ If can be denested itsini* square roots only. Furthermore , w/ien it succeeds y/f will be 
expressed as y/j'= dy/s, for some d e Q h and positive s e Q . 


1 


Proof. The proof is by induction on It. 

Basis (h = 0): ^ has depth 1 over Q, so y/f can be denested if and only il f—d 2 for 
some d e Q . (In this case, s= 1 in the statement of the theorem.) 

Induction (h> 0): Let /= a + by/r h - l9 where a, b e Q h _, and r h _ t e Q. Suppose yZ/can 
be denested over Q using square roots only. Then certainly y/f can be denested over Q h . i 
using square roots only, that is, 


y/a + by/r h - v eQ h - x (y/r k - U y/a x ,... 9y /a k \ 
for some a l9 a 2l . . a k e Q h _ Take K = Q h - { in Theorem 1, and recall that y/r h _ } 
tQ h - x . Hence, applying Theorem 1 , y/a 2 -b 2 r h ~ i e Q h . t . S ince y/a 2 -b 2 r h l has depth 
over Q one greater than any element in Q,,_,, y /a 2 — b 2 r h - l denests over Q using square 
roots only. By the induction hypothesis, the recursive invocation in step (3) will denest 
this correctly as y/a 2 -b 2 r h _ l = ej~s, for some eeQ h - { and positive seQ. 

Since e^fs = y/a 2 -b 2 r h . { eQ h ~i and e e Q h _ ]t y/se Q h - { . The procedure of section 3 
will succeed in expressing y/s as an element of (?/,_!• 

The correctness of 

/ , L / I + 2by/r h - l 

in step (5) is as in Theorem 4. If t e Q, then the right-hand side has depth 1 over Q , so y/f 
has been denested over Q. In this case, the only point to verify is that t is positive. But this 
follows from the fact that 

r t + 2bjr h ., 


2y/a + b y /r h - l 

is real. 

Assume, then, that t e Q h _ t - Q, that is, has depth 2 over Q , so the denesting is not 
completed. Note that _ 

v / t_ f+2fr V^-^ 


Since by hypothesis y/a+bjr h _ \ can be denested over Q using square roots only, and 
since f + 2/) x /r^_ 1 e Q h has depth 1 over Q, y/f can be denested over Q using square roots 
only. By the induction hypothesis, the recursive invocation of step (6) will return j e Q h ^i 
and positive u e Q satisfying y/t—jy/u. The correctness of step (7) is then simple to 
verify. □ 

Consider the example of denesting 

+ hJTt = V(112 + 70^2) + (46 + 34^2)75. 

Here h = 2, a = 112 + 70^2, fc = 46 + 34 N / / 2, and r, = 5. Furthermore, Q 0 = Q, 
Qi = Q(\/ 2), and Q 2 = (?(/, ^/5). The first task is to denest 

V^-^r, = y/(U2 + 70./2) 2 -(46 + 34^/2) 2 x 5 = ^204 + 40^/2 
recursively. This proceeds as in section 2, with the resulting denesting 

d = ^204 + 40^2 = 2+ 10 v /2. 



Step (4) would express tl d-> an element of (2, = (2(,y2), which it already is. In step (5). 
r = 228+160^/2 and 

(228 + 1 60y/l) +2( 46 + 34 v /2) x /5 
2^/228+160^/2 

The numerator is now denested, but J~t in the denominator is not. However, yft contains 
no occurrence of Js, so we proceed to denest J~l recursively in step (6). This proceeds 
again as in section 2, with the resulting denesting 

yft = ^228 + 160^2 = 8 ^/ 2 + 10 . 

Finally, in step (7), 


= (5 + 4 ^ 2 ) + (3 + 72 ) 75 . 

Why does this approach not work if the input formula has depth n > 2 over Q? The 
problem is that ^/t computed in step (5) in general has depth it-1 over Q h - X , which in 
turn will have elements of depth n — 1 over Q. Hence, >w /7 will have depth 2n — 2 over Q , so 
the formula’s depth over Q actually increases if n > 2. 

The authors are indepted to Clifford Earle and Oscar Rothans for the proof of Lemma 5. 



I'yjr~i = 


(228+ 160^/2)+ 2(46+34^/2)75 
2 ( 8 ^ 2 + 10 ) 
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J Simplification of Expressions Involving Radicals 

}' 
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I Algebraically dependent expressions arise in a large variety of symbolic computations. People 

’ seem to have the best intuition about expressions involving radicals. Symbolic computations 

■ with simple, non-nested, radicals is relatively straightforward; however, when the radicals are 
nested the problem becomes more difficult. This paper presents an algorithm for determining a 
linearly independent basis for a set of radicals (nested or not). This allows elementary 
techniques to be used for arithmetic operations on expressions involving elements of this set. 

In addition we provide a structure theorem that provides a sufficient condition for a nested 
radical to be expressed in terms of radicals of lower nesting level. These two techniques are 
powerful tools for computations involving radicals. 

i 

Algebraic manipulation is the study of techniques for performing mathematical 

■ computations symbolically. The most carefully studied and best understood problems in 
' algebraic manipulation have been those that involve polynomials and rational 

functions— expressions involving independent variables combined by addition, 

1 multiplication and division. More complex expressions can often be treated as 
polynomials if the algebraically independent kernels can be identified. For instance, 
e 2x +e x can be treated as a polynomial in the kernel e x . 

; It would be erroneous to view e 2x + e x as a polynomial in e 2x and e x since they are not 
independent. They satisfy the polynomial relation A-B 2 . Expressions that satisfy 
polynomial relations are algebraically dependent. Often more than one kernel is needed to 
■* express an expression as a polynomial over some simple ring. In this case the set of 
, Vemels used is called a basis. Ideally, the elements of the basis should be algebraically 
independent. 

When algebraic numbers (^7) or algebraic functions (y/x 2 -3) are involved things 
, become a little more complex. Algebraic expressions cannot be used as kernels of 
polynomials since they must be reduced when the kernels appear to high enough powers 
if (^-'2 2 -> 2). We are looking for independence results that ensure this is the only 
'j simplification wc need to look for. 

j Regardless of whether the kernels are transcendental or algebraic, there are two basic 
| problems that must be resolved when introducing new types of expressions to an algebraic 
| manipulation system: 

> —Determining a “sufficiently independent” basis for the computation, so that the 
expressions can be treated as polynomials. 

— Determining simple expressions for the basis elements so the user can understand the 
resulting expressions. 
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